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Recently the first chemical systems which show the amplification of enantiomeric excess (ee) was 
found. Inspired by these experiments, we propose a few chemical reaction models in a closed sys- 
tem. The reactions consist of autocatalytic production of chiral molecules and their decomposition 
processes. By means of flow diagrams in the phase space of enatiometric concentrations, chiral 
symmetry breaking is studied. 

Without any destruction processes, the quadratic autocatalysis leads to the ampliflcation of ee, 
but the flow is attracted to a flxed line, which is only marginally stable. The destruction processes 
alter the attractor structure completely. With a cross inhibition the linear autocatalysis is sufficient 
for the chiral symmetry breaking with a complete homochirality. With a back reaction, the quadratic 
autocatalysis is necessary to bring about the chiral symmetry breaking. The final ee is not complete 
if there is a racemization process. 

PACS numbers: 



I. INTRODUCTION 

Because of the three-dimensional tetragonal arrange- 
ment of four chemical bonds around a carbon atom, or- 
ganic molecules can take two types of stereostructures; a 
right-handed one and a mirror-image left-handed oneii 
The molecules which have either of these two forms are 
called chiral, and the two forms are enantiomers with 
each other. In general the two types are classified in R 
and S forms, but, for sugars and amino acids, D and 
L representations arc common in use. From energetic 
point of view, the two enantiomers arc indistinguishable 
and they should exist with equal probability. In fact, 
in the usual chemical reactions both types of enatiomers 
are produced with equal probabilities, and the product 
is called a racemate or a racemic mixture. 

In contrast, life on earth utilizes only one type 
of enantiomers: levorotatory(L)-amino acids and 
dextrorotatory(D)-sugars. This symmetry breaking in 
chirality is called the homochirality, and the origin of 
this unique chirality selection has long been intriguing 
many scientists i Various mechanisms have been pro- 
posed to induce an asymmetry in the primordial molecu- 
lar environment )2iS*^iS*SiL& but the induced chiral asym- 
metry in all the proposed models is expected to be very 
minute. Therefore, the mechanisms to amplify enor- 
mously the enatiometric excess (ee) are indispensable. 

Frank is the first to show that an autocatalytic chem- 
ical reaction with an antagonistic process in open sys- 
tems can lead to an amplification of eeiS Many theoret- 
ical models have since been proposed, but they are of- 
ten criticized as lacking any experimental support i Re- 
cently, a first system which shows the spontaneous am- 
plification of ee has been found by Soai's group >ii*i24iii4 
This enhancement takes place in a closed system, and is 
explained by the second-order autocatalytic reactionii4 
Inspired by these experiments^iitiSiiii^ we constructed 
some simple models of autocatalytic chemical reactions 
in a closed and conserved systemj ^^i^^ and studied con- 
ditions when the chiral symmetry breaking is possible. 



Here we extend these models and study them in detail. 

In §2, we summarize the fundamental aspects of au- 
tocatalytic reaction models where two enantiomers com- 
pete against each other in a closed system for the lim- 
ited resourced The ee amplification is found possi- 
ble if the production of chiral molecules contains the 
quadratic autocatalytic process. However, concentra- 
tions of enantiomers are attracted to a fixed line, which 
is only marginally stable. Therefore, the final config- 
uration might be susceptible to external perturbation. 
In §3, with the mutual antagonism or the cross inhibi- 
tion, the complete homochirality is concluded if there is 
a linear autocatalysis. In §4, with the back reaction, the 
quadratic autocatalysis is found necessary for the chiral- 
ity selection. Summary and discussions are given in the 
last section §5. 



II. AUTOCATALYTIC REACTIONS IN A 
CLOSED SYSTEM 

We start our study about reaction models with au- 
tocatalytic productions of chiral molecules from achiral 
materials. FrankS has considered the chiral symmetry 
breaking in an open chemical system. However, since the 
experiment is performed in a closed system^ we con- 
sider reaction models in a closed and conserved system, 
and explore effects that the autocatalysis alone can lead. 

Suppose two achiral molecules A and B react to pro- 
duce a chiral molecule C. For the purpose of simplicity 
we denote (-R)-C molecule as R and (S')-C molecule as 
S, hereafter. The reactions A-fB-^ R and A-l-B— * S pro- 
ceed in a closed system, and thus they stop when the 
materials A and/or B are exhausted. Here we assume 
that there is an ample amount of material B, and its 
concentration can be regarded as constant all through 
the reaction. Then, the process depends only on the con- 
centrations of A, R and S molecules, which we denote 
a, r and s, respectively. The conservation of the total 
number of participating molecules gives the relation that 
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the total concentration stays constant determined by the 
initial values as 



a + r + s = ao + ro + So = Co = const. 



(1) 



Then the reaction developement can be depicted as a flow 
in the r-s phase diagram. The conservation law (1) limits 
the phase space in the triangular region 



0<r, s<co, < r + s < cq. 



(2) 



As a first step we consider the simplest autocatalytic 
reaction model such that the production from the mate- 
rial A to R or S is described by the following rate equa- 
tions: 



dv 

— = f{r)a = f{r){ca -r-s) 
ds 

— = f{s)a = /(s)(co - r- s) 



(3) 



with 



f{x) = fco + kix + k2X^. 



(4) 

The coefficient fco represents the rate of a nonauto- 
catalitic reaction, ki of a linearly autocatalytic and k2 
of a quadratically autocatalytic reactions. Since all the 
rates fco,i,2 are non-negative, /(x) > for x = r and s. In 
the Soai reaction, the quadratic autocatalysis is shown to 
explain the time evolution of the reaction yield and the 
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A. Fixed points 

We now study the flow diagram in r — s space in typical 
situations. First we search fixed points of Eq.Q in the 
region It is easy to find that a = gives a line of fixed 
points at r + s = Co, which are stable and correspond to 
the state where all the material A is consummed up. Only 
for the case without nonautocatalysis, fco = 0, another 
fixed point is possible at the origin r = ,s = 0, but it is 
easily shown to be unstable. 




FIG. 1: Flow diagrams for (a) nonautocatalytic (fco > 0, fci = 
k2 =0), (b) linearly autocatalytic (fci > 0, fco = ^2 = 0), 
and (c) quadratically autocatalytic (^2 > 0, fco = fci = 0) 
reactions. 

The flow will then be attracted to the fixed line r + s — 
Co, and is given by the equation 

dr f ds 



f{r) J f{s) 



— constant. 



(5) 



Integration of this equation for general f{x) can be car- 
ried out analytically, but we discuss the flow lines only 
in three typical cases in the following. 

• Nonautocatalytic case: 

Without autocatalysis, fco > but ki — k2 = 0, flow 
lines are easily obtained by interation of eq.(5) and they 
have a unit slope, r — s = tq — sq =constant, as shown in 
FigEtta). 

• Linearly autocatalitic case: 

With only a hnear autocatalysis, ki > 0, fco = ^2 = 0, 
the integration of eq.(5) is easy to yield the trajectory 
r/s = tq/sq ^constant. The flow is represented by lines 
projecting out from the origin, as shown in Fig^b). 

• Quadratically autocatalytic case: 

For a quadratically autocatalytic case, fc2 > 0, fco = 
ki = 0, flow lines are a collection of hyperbolae projecting 
out from the origin, 1 / r — 1 / s = 1 / vq — I/so =constant, 
as shown in Fig^c). 



B. Enatiometric Excess 

In order to quantify the ee or the degree of chiral sym- 
metry breaking, we use an order parameter defined as 



(6) 



• Nonautocatalytic case: 

For the non- autocatalytic case in Fig. 1(a), the nu- 
merator r — s remains constant during the time evo- 
lution whereas the denominator r + s increases. Thus 
the ee decreases monotonously from the initial value 
|0o| = K'^o — so)/{ro + so)\ to the asymptotic one l^ool ~ 
{{fo ~ •so)/co|- If the initial concentrations vq and sq are 
small, ee |(/)| becomes very small. This result is quite nat- 
ural from the microscopic point of view, since the nonau- 
tocatalytic reaction produces both enatiomers randomly 
and leads to the racemization. 

• Linearly autocatalytic case: 

For the linearly autocatalytic case in Fig. 1 (b) , the ratio 
r/s remains constant and so does the ee. 

• Quadratically autocatalytic case: 

Only with a nonlinear autocatalysis in Fig. 1(c), the ee 
increases up to some saturation value |0oo|, as shown in 
Fig. 2. This asymptotic value of (jjac depends on the initial 
concentrations tq and sq, even when the initial (f>o and 
the total Co are fixed, as shown in Fig. 2. The smaller the 
initial amount of chiral molecules (tq + sq) are, the more 
the ee is amplified. 

• General case: 

For more general cases, the ee amplification is ex- 
plained by studying the evolution of 0, which is derived 
from eq.(3) as 

§ = ^0h2fco + ^{co - af{l 0^)]. (7) 
dt Co — a 2 

The evolution equation turns out to be independent of 
the linear autocatalytic rate ki. Without quadratic au- 
tocatalysis (fc2 = 0) the ee decreases as we have already 
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FIG. 2: Time evolution of the enantiomeric excess (j) with 
quadratic autocatalytic reaction ^2 > with the same initial 
value 4>o but the different amount of initial R and S molecules, 
To and So. 



observed, since the concentration of A material, a, and 
the sum of R and S products, co — a = r + s, should 
always be positive. The situation does not alter by in- 
cluding only a linear autocatalysis with fci > 0: it does 
not directly affect the evolution of (f). Only when the 
quadratic auatocatalysis dominates at > 4A:o/co, the 
ee can increase. Therefore, in this simplest model only 
the quadratic autocatalysis can enhance the chiral assym- 
metry. This may correspond to the situation observed in 
the Soai reaction. It has to be noted that the asymp- 
totic value (j)oo cannot be determined from eq.(7), since 
lim a(t) = in the present model. 

t — *oo 

Even though the amplification of the enatiometric ex- 
cess or the amplification of the chiral assymmetry can be 
brought about by the nonlinear autocatalysis, the final 
state of the system is on the line of fixed points, and 
thus the final values of ee can be anything depending on 
the initial condition. Since these points are marginally 
stable, i.e. they have no quaranteed stability along the 
direction of the fixed line, they may be susceptible to ad- 
ditional perturbations. We now consider the effect of two 
types of destructive reactions of chiral molecules. 



A. Fixed points 



Fixed points have to satisfy the relation 



f{r)a = /(s)a = /irs. 



(9) 



The simplest solution is a = and rs ~ 0, which cor- 
responds to the end corners of the triangular region, 
(r, s) = (co, 0) and (0, Co). These are easily shown to 
be stable with complete homochirality, = ±1. Other 
fixed points are also possible with a nonzero a and 
/(?■) = /(s) = jirs/a. We shall study them in various 
cases below. 

• Nonautocatalytic case: 

Only with a nonautocatalytic reaction (fcg > 0, fci = 
^2 = 0), f{r) = /(s) = /cq E^nd there is a fixed line 

(,+ ^)(, + ^) = ^(co + ^) (10) 

passing through the two terminal points (r, s) = (cq, 0) 
and (0, Co). This line is an accumulation of infinitely 
many marginal fixed points. In fact the fiow trajectory 
is easily integrated as r — s = /'o — sq ^constant, the same 
lines with a unit slope obtained in the previous section 
without cross inhibition. The cross inhibition only shifts 
the fixed line from the border r -I- s = cq in Fig. 1(a) to 
the hyperbola (tlT 
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FIG. 3: Flow diagrams with cross inhibition /icq = 1 at (a) 
fco > 0, kx=k2=Q and (b) fci > 0, fco = ^2 = 0. 



III. WITH CROSS INHIBITION 

Frank has introduced the mutual antagonism between 
the enantiomers, R and S, and succeeded in explaining 
the chiral symmetry breaking in an open reaction system. 
Here we extend his model to the closed system. The reac- 
tion equation now incorporates terms of cross inhibition 
as 

^ = f{r)a - firs 
ds 

— = f{s)a-- firs. (8) 



• Linearly autocatalytic case: 

With a finite fci as well as fco, a fixed line is replaced 
by a symmetric fixed point U at 



-2fco + fcico + J {2ko + fciCo)2 -I- 4koCofi 

r = s = ^ . (11) 

2(/i + 2fci) ^ ' 

At fco = 0, there appears an additional fixed point at 
the origin r = s = 0. All these fixed points are easily 
shown to be unstable, as shown in Fig.3(b). The only 
stable fixed points are those at the corners, (co,0) and 
(0, Co), with the complete homochirality. Thus, with only 
a linear autocatalysis and without any nonlinear ones, 
complete homochirality is possible with a cross inhibition 
even in a closed system. 

With quadratic autocatalysis, the overall structure of 
the flow is similar to the case with a linear autocatalysis. 
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B. Enantiomeric excess 

Up to the linear autocatalysis (fc2 = 0) with a cross 
inhibition (/i 7^ 0), the ce foUows the evolution equation 

^ = 0[-2fcoa+g(co-a)2(l-(/)2)]. (12) 

at cq — a L 

The rate fci of the linear autocatalytic process is not in- 
volved in the symmetry breaking effect of (\) neither in the 
present case with cross inhibition nor in the case with- 
out cross inhibition in the previous section. One notices 
the similarity between the evolution equation H12() and 
eq.Q derived in the previous section. The effect of the 
quadratic autocatalytic reaction, hi , in eq. Q) is replaced 
by ^/a in ea. H12() . In the present case, the symmetry re- 
covering effect fco weakens as the material A is consumed 
(a 0), and the complete homochirality </) = ±1 is 
achieved asymptotically. For example, the time evolution 
of r/co, s/cq and 4> is depicted in Fig. 4 at fcg = 0, fciCo = 
1,/iiCo = 1 and ro/cQ = 0.0714, sq/co = 0.0700, (/)o = 0.01. 
Initially r and s increase simultaneously to approach the 
symmetric fixed point U: (r/cg, s/co)=(l/3,l/3) with a 
constant <j) ~ 4>q, driven by the linear autocatalytic pro- 
cess ki. However, since U is unstable, the trajectory is 
eventually expelled from U and r increases at the cost of 
s. As a result (j) increases to the full homochirality. 
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FIG. 4: Time evolution of r,s, and cj) ^t ko — 0, fcico 
1, ^co = 1 and ro/co = 0.0714, sq/cq = 0.0700, 0o = 0.01. 



IV. WITH BACK REACTIONS 

We now consider another type of decomposition pro- 
cess, a simple back reaction from the product R or S to 
an A molecule. Governing equations are 



f{r)a - Xr 
f{s)a - As 



(13) 



Since the back reaction is linear in the product concen- 
trations r and s, the evolution of the ee remains the same 
with eq.(7). Therefore, the quadratic autocatalysis with 
a finite ki is the necessity to break the chiral symmetry in 
the present case. Even though follows the same evolu- 
tion (7) without the back reaction, its asmptotic behavior 



is quite different from that in §2, because the fixed point 
structure is drastically altered by the back reaction. 
Fixed points 

Fixed points of the flow satisfy conditions 



f{r)a - Xr = f{s)a - Xs = 0. 



(14) 



It is trivial that a = cannot be a fixed point any more 
unless Co = 0. As we have supposed, the fixed line r + s = 
Co in §2 is very sensible to the perturbation like a back 
reaction. We now study structures of fixed points in some 
simple cases. 

• Nonautocatalytic case: 

In the non-autocatalytic case with fco > 0, fci = fc2 = 
with a finite A, there is only a symmetric fixed point at 



Co/2 



< 



Co 



(15) 



1 + A/2fco ~ 2 
which is easily shown to be stable, as shown in Fig. 5(a) 




FIG. 5: Flow diagrams for (a) fco > 0, fci = fc2 = 0, A = 1, (b) 
fci = 1, fco = fc2 = 0, A = 0.1, (c) fc2 = 1, fco = fci = 0, A = 0.1, 
(d) fe2 = 1, fco = fc2 = 0,A = 0.2, (e) fcico = 0.4, fcacg = 
1, fco = 0, A = 0.3. (f) fco = 0.05, fciCo = 0., fc2cg = 1, A = 0.1. 




FIG. 6: Time evolution of the ee (p with nonlinear autocatal- 
ysis and a back reaction (a) without random creation fco = 0, 
and (b) with a finte fco = 0.05. 

• Linearly autocatalytic case: 

With only a linear autocatalysis fci > 0, fco = ^"2 = 0, 
we have only a trivial fixed point at the origin (r, s) = 
(0,0), when the rate of back reaction is as large as 
A > fci Co: there will be no net production of chiral 
C molecules because of strong back reaction. On the 
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contrary, when the rate of back reaction is small as 
A < fciCo, there are additional fixed points on a line 
r + s = Co — A/fci. This line corresponds to the one 
r + s = Co in §2 shifted by the back reaction A. In fact 
the trajectory is the same as in §2 and radiates from the 
origin r/s = tq/so ^constant, as shown in Fig. 5(b). 
• Quadratically autocatalytic case: 

With a quadratic autocatalysis fc2 > 0, fco = fci = 0, 
the back reaction decimates the fixed line r + s = co and 
produces several fixed points. One of them is a trivial 
fixed point at the origin (r, s) = (0,0), which is of no 
interest. 

(1) With a small rate of back reaction, A < fc2Co/8, 
there are other six fixed points: two on the r axis, 
Ui : (A_,0) and Si : (A+,0), two on the s axis, 
U2 : (0,^-) and S2 : (0,^+), and two symmetric 
ones, U3 : {B^,B^) and U4 : {B+,B+), with 

As an illustration shown in Fig. 5(c), we use pa- 
rameter values A:2C§ ~ 1, A = 0.1, which gives 
fixed points Ui_4 and Si. 2 specified by the values, 
A+ = 0.8873CO = Co - A_, and B+ = 0.3618co = 
Co/2 — B-. The stable fixed points Si and S2 cor- 
respond to the states with complete homochirality 
with = ±1. Except the region around the trivial 
fixed point O, the phase flow is attracted to one of 
these two homochiral states, Si or S2. Generally, 
the region of initial points which will eventually be 
attracted to a certain fixed point is called the basin 
of attraction. In the present case, the r — s phase 
space is divided into three basins of attractions to 
Si, S2, and O, by separatrixes U3U1, U3U2, and 
U3M. 

The time evolution of the order parameter (p is de- 
picted in Fig. 6 (a) with the parameters kg = ki = 
0, A;2Co = 1, Aco = 0.1, the same ones used to obtain 
the flow diagram in Fig. 5(c). The initial condition 
is To = 0.2 + 0.001, So = 0.2. The flow is almost 
symmetric as r w s first and approaches to the 
symmetric fixed point U4, but since U4 is unstable 
r starts to increase at the cost of s and eventually 
(j) approaches the full value, unity. 

(2) As the back reaction rate increases, two symmetric 
fixed points U3 and U4 merge to disappear for A > 
^2Co/4. There remains two homochiral fixed points 
Si and S2. Their basins of attraction are close to 
the r or s axes. Close to the diagonal r = s the flow 
runs to the trivial fixed point O. The separatrixes 
of three basins of attractions are UiMi and U2M2. 
With the parameter values fc2Co = 1, A = 0.2, fixed 



points Si and S2 are determined with parameters 
= 0.7236co and A- = cq — A^, as shown in 
Fig.5(d). 

(3) Finally, when the back reaction is too strong as 
A > fc2Co/4, all the flow runs to the origin. 

• Linear and quadratic autocatalysis: 

We consider now the effect of a linear autocatalytic 
process in addition to the quadratic autocatalysis ^2- As 
far as ki is smaller than A/co, the origin r ~ s = 
is attractive, and the schematics of the flow in the r — s 
phase space is the same with those for fci = 0. When ki is 
larger than A/cq, the origin r = s = looses the stability, 
and there are only three fixed points: two stable ones on 
r and s axes Si: (C, 0) and S2: (0, C) and a symmetric 
and unstable one at U: {D, D) with 

_ fc2Co - ki + ^/ (fc2Co + fci)^ - 4fcjA 
^ 2fc2 

^ _ fc2Co - 2fci + A/(fc2Co 4- 2fci)2 - 8fc2A 
4/C2 

and the typical flow in this case with ko = 0, fcico = 
0.4, fc2Co = 1, A = 0.3 is shown in Fig. 5(e), which looks 
similar to the case with the cross inhibition shown in 
Fig. 3(b): complete homochiral states, (j) = ±1, is at- 
tained asymptotically even though not all the A material 
is consumed. 

• Autocatalysis with nonautocatalytic production 

We found that the quadratic autocatalysis with back 
reaction permits to attain the homochiral situation, if 
there are a few C molecules initially outside the basin of 
attractive to the origin r = s = 0. In order to actuate 
the autocatalysis, however, one has to provide initial C 
molecules. Therefore, nonautocatalytic production with 
a finite kg seems necessary. We now consider cases when 
the random creation of racemic mixtures by fco coexist 
with the autocatalytic production as well as the back 
reaction. 

(1) With a linear autocatalysis fco > 0, fci > 0,fc2 = 
0,A > 0, the fixed line for fco = in Fig. 5(b) is 
replaced by a symmetric fixed point at 

r = s =-^(fciCo - 2fco - A 
4fci L 

+ V(fcico - 2fco - A)2 + 8fc-ofc-ico)} (16) 

which is stable and = 0. Random racemization 
effect by fco dominates and recovers the symmetric 
state. 

(2) With a quadratic autocatalysis fco > 0, fc2 > 
0, fci ~ 0, the situation is rather complex. Since 
/(O) ^ 0, all the fixed points stay within the trian- 
gle (2); a, r, s ^ 0. Then the fixed points satisfy 



6 



the relation 

ko + k^r"^ _ A _ fco + ^25^ 
r a s 

{r ^ s){kQ ^ k2rs) = (17) 

The symmetry broken states stay on a hyperbola 
rs ~ ko/k2. They satisfy further the relation 



Aa(r) = k 



1° 
4 

^ A > 0. 



Co 

2 



fcp 
k2r 



(18) 



If fco is smaller than fc2Co/16, the function Aa{r) in 
eq.(18) has two peaks with a height Xmax ~ fc2Co/4, 
as shown in Fig. 7(a). It also has a local minimum 

value Amin = 2y/%{coy/k^ - 2y/%) at r = A/fco/fc2- 

Then for small A < Amm there are two solutions of 
asymmetric fixed points which are stable, for inter- 
mediate case with Xnmi < A < Xmax there are four 
solutions two among them are stable and two other 
are unstable, and for large A > X^ax there is none. 

If fco lies between fc2CQ/16 and fc2C§/4, Aa(r) 
has only a single peak with a value AJ,j^^ = 
2\/k(){caVk2 — 2^/ko). Therefore the number of so- 
lutions are two for small A < A^^^ which are stable, 
and none for large A. 

Further larger fco > fc2Co/4, there is no asymmetric 
fixed point for any A. The last result is already 
drawn by the time evolution of </> given in Eq.(7), 
which is valid in the present case, too. 
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FIG. 7: Functions (a) Aa{r) and (b) As(r) to determine (a) 
asymmetric and (b) symmetric fixed points. 

As for the symmetric solutions with r ~ s, the 
effect of random creation fco is more effective. Since 
the fixed points have to satisfy a relation 



A.(r) 



i-2Co 



1\3 



if---) 



Co 

fco ]_\ 

2cg 108/ 



^2 eg 

12 



(19) 



the number of fixed points reduces to one which is 
unstable if fco is larger than fc2CQ/108, as shown in 



Fig. 7(b). For a small rate fcg < fc2CQ/108, the num- 
bers of symmetric fixed points depends on A. For a 
small A, there is one unstable, for an intermediate 
value of A there are three among which one is stable 
and two unstable, and for a large A there remains 
only one stable one. For example, by adding fcg 
as small as fco = 0.05(> fc2CQ/108) to the situation 
shown in Fig. 5(c), the origin looses the stability and 
there remains only a single symmetric fixed point 
which is unstable. Therefore, the flow diagram re- 
duces to the one shown in Fig. 5(f). The achieved 
ee is not complete, but is very large if fco is small, 
as shown in Fig. 6(b). There the initial condition 
is chosen as ro = 0.1 + 0.001 and s = 0.1, but the 
final (f> deos not depend on the initial values of ro 
or So- 



V. SUMMARY AND DISCUSSIONS 

We studied various autocatalytic models with back re- 
action or with cross inhibition processes to search for the 
conditions to achieve the chirality selection in a closed 
system. The results arc summarized in the Table 1. 
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r ~ (fco + kir + k2r^)a — Ar — /irs 
s = (fco + kis + k2S^)a — As — firs 
a — Co — r — s 



TABLE I: Classification of reaction models witli parameters 
defined in the bottom column. Symbol > means that the 
value is non-negative and > means that the value must be 
positive. B means that the final state is chiral symmetry 
breaking and S means otherwise. Details for each entry are 
explained in the indicated sections. 



If the production of chiral molecules R and S from 
the achiral material A involves a linear autocatalysis, the 
cross inhibition between R and S is necessary to realise 
the homo chirality, as is summarized for the cases (l)-(4) 
in Table 1. The case (3) corresponds to an extension of 
Frank's model to a closed chemical system. 

With a quadratic autocatalysis, the amplification of 
enantiomeric excess (ee) is achieved without any destruc- 
tion process, but the steady state can have any finite 
value of ee ( — 1 < (p < 1), as the stable fixed line exists 
extending from cj) = 1 to (j) = —I. In fact the final value 
of ee depends on the details of the initial values of chiral 
as well as achiral materials' concentrations. Therefore 
we do not regard this case (1) as the chiral symmetry 
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breaking in a restricted sense. 

With a linear back reaction process and a quadratic 
autocatalysis as a special case of (5), the complete ho- 
mochirality is possible to realize when the rate of back 
reaction is sufficiently small. The minority enantiomer 
is decomposed back to the material A and recycled to 
produce the majority enantiomer and vice versa. This 
decomposition process works more effectively for the mi- 
nority because of the suppressed production due to the 
nonlinear autoctalysis. 

However, with only autocatalytic processes, the initial 
production of chiral molecules is impossible. Therefore, 
the random production of chiral molecules by the nonau- 
tocatalytic process described by a nonzero value of coeffi- 
cent fco is indispensable. For a small value of fco, the phase 
space structure of the flow and fixed points is not qual- 
itatively affected and resembles with that of vanishing 
fco case. However, for fco greater than the critical value 
fc2Co/108, this process affects the structure of racemic 
fixed points so drastically that the system can have only 
almost homochiral states as stable fixed points, as shown 
in Fig.5(f). 

We have studied two types of simple models (3) and 
(5) where the chiral symmetry can be broken in the 
closed system and complete or almost complete homochi- 
ral states can be achieved. One type of the model with 
the cross inhibition (3) seems powerful for this purpose, 
but the chemical system which shows this process has not 
been found so far and it is very interesting to find actual 



chemical system. On the other hand, the nonlinear auto- 
catalytic system is found by Soai's group and thus, there 
is a real possibility for a nonlinear autocatalytic model 
with back reaction (5) to exist. But we are not certain 
whether the back reaction proceeds in a measurerable 
time for the experiment. 

From these analysis of present simple models based 
on closed systems, we would like to venture to specu- 
late a possible scenario for the initial preparation of con- 
stituent materials of life as follows. In a confined nook 
and cranny somewhere in the premordial sea, there was a 
very favorable conditions for the autocatalytic reactions 
to take place and after a sufficiently long time, the initial 
homochiral selection was established. Once this selec- 
tion was completed in this tiny closed corner, the chiral 
molecules began to spread out and to conquer eventually 
the surrounding region and futher all over the earth. 

So far, we have assumed that the system is homoge- 
nious, such that the whole system is described by only 
a space-independent concentrations. In the actual situa- 
tions, the concentrations should vary in space. We have 
introduced a lattice model for the reaction-diffusion sys- 
tem and studied the spreading of chiral symmetry break- 
ing. The interested readers are invited to refer to the 
original paper ^ 
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